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ABSTRACT 

Models with late time cosmic acceleration, such as the A-dominated CDM model, 
predict a freeze out for the growth of linear gravitational potential at moderate redshift 
z < 1, what can be observed as temperature anisotropies in the CMB: the so called 
integrated Sachs- Wolfe (ISW) effect. We present a direct measurement of the ISW 
effect based on the angular cross-correlation function, wtg(9), of CMB temperature 
anisotropies and dark-matter fluctuations traced by galaxies. We cross-correlate the 
first-year WMAP data in combination with the APM Galaxy survey. On the largest 
scales, 8 = 4—10 deg, we detect a non- vanishing cross-correlation at 98.8% significance 
level, with a 1-a error of wtg — 0.35 ± 0.14//K, what favors large values of fl\ ~ 
0.8 for flat FRW models. On smaller scales, 6 < fdeg, the correlations disappear. 
This is contrary to what would be expected from the ISW effect, but the absence 
of correlations may be simply explained if the ISW signal was being cancelled by 
anti-correlations arising from the thermal Sunyaev-Zeldovich (SZ) effect. 



1 INTRODUCTION 



The recent measurements of CMB anisotropies made public 
by the WMAP team are in good agreement with a 'concor- 
dance' cosmology based on the ACDM model. The unprece- 
dented sensitivity, frequency and sky coverage of this new 
data set provides us the opportunity of asking new questions 
about the evolution of the universe. In this paper we present 
the cross-correlation of the cosmic microwave background 
(CMB) anisotropies measured by WMAP (Bennett et al 
2003a), with galaxies in the APM Galaxy Survey (Maddox 
et al 1990). In the observationally favored ACDM mod- 
els, a non-vanishing CMB-galaxy cross-correlation signal 
arises from the distortion of the pattern of primary CMB 
anisotropies by the large-scale structures as microwave pho- 
tons travel from the last scattering surface to us. On large 
angular scales such distortion is mainly produced by the 
energy injection photons experience as they cross time- 
evolving dark-matter gravitational potential wells: the so- 
called integrated Sachs-Wolfe effect (Sachs & Wolfe 1967). 
On smaller scales additional secondary anisotropies are pro- 
duced when photons scatter off hot intra-cluster gas, ie the 
Sunyaev-Zeldovich effect (Sunyaev & Zeldovich 1969). 

In this work, we shall use optical galaxies from the 
APM survey as tracers of the large-scale dark-matter dis- 
tribution of the universe. The APM Survey has produced 



one of the best estimates of the angular galaxy 2-point cor- 
relation function to date. Its shape on large scales led to 
the discovery of excess large-scale power, and gave early in- 
dications of the ACDM model (Efstathiou, Sutherland, & 
Maddox 1990; Maddox et al 1990; Baugh & Efstathiou 
1993), (Gaztanaga 1995). Higher-order correlations have 
also been studied in the APM Galaxy Survey (Gaztanaga 
1994; Szapudi et al 1995; Frieman & Gaztanaga 1999). For 
the first time, these measurements were accurate enough 
and extended to sufficiently large scales to probe the weakly 
non-linear regime with a reliable Survey. The results are in 
good agreement with gravitational growth for a model with 
initial Gaussian fluctuations. They also indicate that the 
APM galaxies are relatively unbiased tracers of the mass on 
large scales (Gaztanaga & Juszkiewicz 2003). Moreover the 
APM results are in excellent agreement with other wide field 
photometric surveys, such as the Sloan Digital Sky Survey 
(SDSS), for both number counts and clustering (Dodelson et 
al. 2002; Scranton et al 2002; Gaztanaga 2002a; Gaztanaga 
2002b). 
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2 DATA & SIMULATIONS 

The APM Galaxy Survey (Maddox et al 1990) is based 
on 185 UK IIIA-J Schmidt photographic plates each corre- 
sponding to 5.8 x 5.8 deg 2 on the sky limited to bj ~ 20.5 
and having a mean depth of ~ 400 Mpc/h for b < —40 
deg and 8 < — 20 deg. These fields were scanned by the 
APM machine and carefully matched using the 5.8 x 0.8 
deg 2 plate overlaps. Out of the APM Survey we considered 
a 17 < bj < 20 magnitude slice, which includes 1.2 million 
galaxies at a mean redshift z = 0.15, in an equal- area pro- 
jection pixel map with a resolution of 3.5', that covers over 
4300 deg 2 around the SCO 

We use the full-sky CMB maps from the first-year 
WMAP data (Bennett et al 2003a). In particular, we have 
chosen the V-band (~ 61 GHz) for our analysis since it 
has a lower pixel noise than the highest frequency W-band 
(~ 94 GHz), while it has sufficiently high spatial resolution 
(21') to map the typical Abell cluster radius at the mean 
APM depth. We mask out pixels using KpO mask, which 
cuts 21.4% of sky pixels (Bennett et al 2003b), to make sure 
Galactic emission does not affect our analysis. WMAP and 
APM data are digitized into 7' pixels using the HEALPix 
tessellation 1 . Figs show these maps smoothed using a 
Gaussian beam of FWHM = 5 deg (left) and 0.7 deg (right 
panels) . 

To determine the accuracy of our error estimation we 
run 200 WMAP V-band Monte-Carlo realizations. We sim- 
ulate the signal by making random realizations of the CMB 
angular power-spectrum as measured by WMAP, convolved 
with its measured symmetrized beam profile, to which we 
add random realizations of the white noise estimated for 
the V-band (Hinshaw et al 2003). Sampling variance in the 
WMAP-APM cross-correlation is thus evaluated by com- 
puting the correlation between the simulated V-band CMB 
maps (with WMAP KpO mask pixels removed) and the APM 
survey. 



3 WMAP-APM CROSS-CORRELATION 

We define the cross-correlation function as the expectation 
value of density fluctuations 5g = No/ < No > — 1 and 
temperature anisotropics At = T — To (in /xK) at two po- 
sitions hi and fi2 in the sky: wtg{9) = (AT(ni)<5G(n 2 )), 
where 9 = |ri2 — fii|, assuming that the distribution is sta- 
tistically isotropic. To estimate wtg{9) from the pixel maps 
we use: 



wtg(8) = 



E^.AT(rii) fcfa) 



(1) 



where the sum extends to all pairs i, j separated by 9 ± A9. 
The mean temperature fluctuation is subtracted (At) = 0. 
The weights Wi can be used to minimize the variance when 
the pixel noise is not uniform, however this introduces larger 



1 Some of the results in this paper have been 
derived using H EALPix (Gorski et al 1998), 
http: / /www. eso.org/science/healpix 



cosmic variance. Here we follow the WMAP team and use 
uniform weights (i.e. Wi = 1). We consider angular scales, 
9 < 10 deg. Cross-correlations are expected to be dominated 
by sampling variance beyond ~ 10 deg, where the APM an- 
gular correlation function vanishes. Fig|5|shows the resulting 
cross-correlation. On scales above 9 > 5 deg there is a sig- 
nificant correlation above the estimated error-bars. 

Fig HI shows the 1-a confidence interval for wtg{9) ob- 
tained using the jack-knife covariance matrix. Surveys are 
first divided into M = 10 (similar results are found for 
M = 10 — 20) separate regions on the sky, each of equal 
area. The analysis is then performed M times, each time 
removing a different region, the so-called jack-knife subsam- 
ples, which we label k = 1 . . . M. The estimated statistical 
covariance for wtg at scales 9i and 9j is then given by: 



M 1 M 
dj = — ^ A^tg(^) Aw tg(^j 



(2) 



where Atu TG (0j) is the difference between the fc-th subsam- 
ple and the mean value for the M subsamples. The case i = j 
gives the error variance. The accuracy of the jack-knife co- 
variance have been tested for both WMAP (Gaztanaga et al 
2003) and the APM and SDSS survey (Zehavi et al. 2002; 
Scranton et al 2002; Gaztanaga 2002a; Gaztanaga 2002b). 

We have used the Monte-Carlo (MC) simulations de- 
scribed above to find that the jack-knife (JK) errors from 
a single MC simulation agree very well (better than 20% 
accuracy), with the MC error from 200 realizations. This 
is shown in Fig. The JK error from a single realization 
(shown as squares) is in excellent agreement with the MC 
error (dispersion from 200 realizations, shown as continu- 
ous line). Thus, in our case, using JK errors over a single 
realization gives an unbiased estimation of the true error, 
but the variance (shown as errorbars in Fig. [3J in this error 
estimation can be as large as 20%. 

On the other hand, the JK errors in the real WMAP- 
APM sample are shown as a dashed line in Fig- 01 They are 
comparable to the errors we measure using the MC simula- 
tions on scales 9 > 4 deg, where the JK errors are potentially 
subject to larger biases (as we approach the size of the JK 
subsample). On smaller scales, the JK errors from WMAP 
are up to a factor 2 smaller than the JK errors (or sample 
to sample dispersion) within the MC simulations. The de- 
viation is significant, given the errorbars, and it therefore 
shows that the MC simulations (rather than the JK error 
method) fail to reproduce the WMAP-APM data on small 
scales. This is not very surprising as the MC simulations do 
not include any physical correlation that might be present 
in WMAP-APM and assume a CMB power spectrum that 
is valid for the whole sky, and not constraint to reproduce 
the WMAP power over the APM region. Alternatively, the 
JK errors provide a model free estimation that is only sub- 
ject to moderate (20%) uncertainty, while MC errors depend 
crucially on the model assumptions used to produce simula- 
tions. 

Despite these differences in the errors on small scales, 
the overall significance for the detection on different scales is 
similar when we use the \ 2 values below or when we ask how 
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Figure 1. APM galaxy density fluctuation maps (top panels) compared to WMAP (V-band) maps (bottom panels) and the cross- 
correlation map (middle). In each case left panels show the maps smoothed with a Gaussian beam of FWHM = 5 deg, while right panels 
have FWHM = 0.7 deg. We use normalized units (dimensionless) with linear color scheme in the range (— 3<r, +3cr), being a the pixel 
variance in each map. 



many of the 200 MC simulations have a signal and JK er- 
ror comparable to the observations. When a particular MC 
simulation has an accidentally large cross-correlation signal, 
it also has a large noise (JK error) associated. We find that 
at 5 degrees only one of MC simulations have a signal to 
noise ratio comparable to the observed WMAP- APM corre- 
lation on large scales. This means that the significance of the 
cross-correlation detection is better than 1%. We can now 
estimate confidence regions in wtg using a x 2 test with the 
JK covariance matrix: 



N 



where Aj 



Ai C- 1 A,, (3) 



ra(@i) ~ w TG(@i) i s the difference between the 
"estimation" E and the model M. The x 2 test gives a best 
fit constant wtg = 0.35 ± 0.13pK (the error corresponds to 
68% C.L. Xmin = 0-19) for the 3 bins in the range 6 = 4-10 
degrees. For this range of scales, a constant null correlation 
wtg = gives Ax 2 =6.1 which has a probability of P = 
1.2% and sets the significance of our detection to 2.5a. We 
find similar results when we use the MC covariance matrix 
(eg, Xmin = 0.33 for best fit; Ax = 5.3 and P = 2% for 
the significance of the detection) but concentrate on the JK 
matrix from now on for the reasons given above. 
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Figure 2. Comparing measurement to predictions: The two solid 
lines and the dotted line show wtg results for WMAP bands 
V, W, and the foreground "cleaned" map. Boxes give the 68% 
confidence levels. Also shown are the theoretical predictions for 
ISW and SZ (upper and lower short-dashed lines) and their sum 
(long-dashed line) for the best-fit model with — 0.8. 
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Figure 3. Errors in the cross-correlation wtg{8) from the dis- 
persion in 200 Monte-Carlo simulations (solid line) as compared 
with the mean and dispersion (squares with errobars) in the jack- 
knife error estimation over the same simulations. Dashed line cor- 
respond to the jack-knife error in the real WMAP-APM sample. 



3.1 Comparison with Predictions 

Galaxy fluctuations in the sky can be modeled as, Sa(n) = 
J dz ()>g{z) 5g{ji,z) where 4>g(z) models the survey selec- 
tion function along the line-of-sight. We will assume here 
that APM galaxies are good tracers of the mass on large 
scales (see so that we can use the linear bias rela- 
tion: 8g — b8, with b ~ 1 and for the power spectrum: 
PGik,z) — b(z) 2 P(k, z). In the linear regime we further 
have: P(k,z) = D 2 (z) P(k). We can then define a galaxy 
window function Wq{z) ~ b(z) D(z) 4>g(z) accounting for 
bias, linear growth, and the galaxy selection function. For 
the APM selection we use the function in (Baugh & Efs- 
tathiou 1993) with a mean redshift z ~ 0.15. Thus the 
galaxy 2-point angular correlation is (Baugh & Efstathiou 
1993),(Gaztafiaga & Baugh 1998) w G a{0) =< 8 G 5 G >= 
J dk k P(k) gc{k9), where the kernel gc(k9) is a line-of- 
sight integral, ga(k9) — l/2w J dz W G (z) jo(k9r) where 
jo is the zero-th order Bessel function, and r(z) denotes co- 
moving distance. We use (Bond & Efstathiou 1984) for the 
linear power spectrum, with shape parameter P = Q m h, and 
h = 0.7. We take as — 1 and b — 1 which give a resonable 
match to the measured variance in the APM (Gaztanaga 
1995). 

The temperature of CMB photons is gravitationally 
blueshifted as they travel through the time-evolving dark- 
matter gravitational potential wells along the line-of-sight, 
from the last scattering surface z„ = 1089 to us today, 
z = (Sachs & Wolfe 1967). At a given sky position n: 
AT ISW (n) = -2 J dz i(n,z). For a flat universe V 2 $ = 
— AnGa 2 p m 8 (Peebles 1980), which in Fourier space reads, 
$(k,z) = -3/2Q m (H /kf S(k,z)/a. Thus: 



ISW to \ 

w TG (9) 



< ^ T SW 8 G >-- 



dk 



P(k) g(k8) 



(4) 



g(k9) = — I dzW IS w(z)W G (z) j (k9r) 



where the ISW window function is given by Wisw = 
-3 Q m (^) 2 F(z), with c/H ~ 3000ft Mpc" 1 the Hub- 
ble radius today. F = d(D/a)/dr = (H/c)D(f - 1) with 
/ ~ Qm 11 (z) quantifies the time evolution of the gravi- 
tational potential. Note that F decreases as a function of 



increasing redshift (as Q m (z) — > 1). It turns out that for 
flat universes, Q m + Qa = 1, Wisw has a maximum at 
Qa — 0.6 and tends to zero both for J1a — » 1 (since the pref- 
actor S} m — > 0) and also for J7a — ► (because F — > 0). The 
ACDM prediction with Qa — 0.8 and T — 0.14 is shown as 
a short-dashed line in Fig 01 

For the thermal Sunyaev-Zeldovich (SZ) effect, we just 
assume that the gas pressure S gas fluctuations are traced 



by the APM galaxy fluctuations S gas 



5 G with an 



amplitude b gas ~ 2, representative of galaxy clusters. Note 
that analytical results based on halo models and hydrody- 
namic simulations show that this "gas bias" factor is scale 
and redshift dependent (Refregier & Teyssier 2002). How- 
ever, for low-2 sources and linear scales one can safely take 
b g as = 2—4. Note that the cross-correlation function is domi- 
nated (its overall shape and amplitude) by large-scale modes 
that are well described by linear theory, but a more precise 
calculation requieres non-linear corrections. Thus a rough 
conservative estimate is given by (Refregier et al 2000): 
Wtc(6) = -b gaa ATw G g(6>), where AT ~ 6.6^K is the 
estimated mean SZ fluctuation from APM clusters, what 
corresponds to a Compton parameter y ~ 2 x 10~ 7 (Re- 
fregier et al 2000). The SZ prediction along with the total 
correlation, wtg = Wtg + w tgi are given by the lower 
short-dashed and long-dashed lines in Fig [5] Deriving more 
accurate parameter contraints from the SZ effect requires 
including non-linear effects in the power spectrum, what is 
beyond the scope of this paper. 



4 DISCUSSION 

The main result of this paper is a measurement of a positive 
cross-correlation wtg = 0.35 ±0.13/^K (l-cr error) between 
WMAP CMB temperature anisotropies and the Galaxy den- 
sity fluctuations in the largest scales of the APM galaxy sur- 
vey, 9 ~ 4 — 10 deg. The significance of this detection is at 
the 98.8% confidence level (2.5 a). Large-scale modes from 
the primary SW temperature anisotropies introduce large 
sampling variance and makes measurements of the ISW con- 
tribution intrinsically noisy. The measured cross-correlation 
on 9 > 4 deg scales is in good agreement with ISW effect pre- 
dictions from standard ACDM models. Using the theoretical 
modeling in §111. A we find a 2-a interval of Qa = 0.53 — 0.86, 
with a best fit value of Qa — 0.8. 

If the detected cross-correlation is only due to the 
ISW effect (Crittenden & Turok 1996), one would expect 
a stronger ISW-induced correlation on smaller scales (see 
FigHJ. Instead, on scales 0.2 < 6 < 1 deg, the mean cross- 
correlation becomes negative, wtg ~ —0.06 ± 0.16/iK. This 
sudden drop can be understood as thermal SZ contribution 
from hot gas in galaxy clusters (Refregier et al 2000). The 
SZ effect contributes to a level w T g = wtg — ii'tc*' ~ 
—0.41 ± 0.13 ± 0.16/xK, where errors reflect the uncertain- 
ties at large and small scales. This result can be used to set 
bounds on the mean Compton scattering of CMB photons 
crossing clusters (see also (Diego, Silk & Sliwa 2003)). De- 
spite the high galactic latitude (b < —40 deg), our results 
can potentially be contaminated by Galactic dust (Nolta 
et al 2003). However, as illustrated in Fig 2, using the 
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KpO masked W-band, V-band or a foreground "cleaned"map 
(Tegmak et al 2003) all give similar results within the errors. 

Fig Q shows the ISW and SZ contributions at the 
map level. On larger scales the APM-WMAP product map 
shows a clear correlation with the APM structures, while on 
smaller scales (right panels) this correlation fades away or 
turns into anti-correlation at the core of the largest struc- 
tures in the APM map. These APM structures correspond 
to the very large scale potentials hosting superclusters or a 
few large clusters in projection. Some of them appear to be 
anti-correlated with WMAP in the product map, but with 
very similar shapes (regions pointed by an arrow in Fig0. 
This is not surprising since the measurement of the ISW ef- 
fect is intrinsically affected by sampling variance from the 
larger amplitude modes due to primordial SW fluctuations. 
In particular, if a real (positive amplitude) ISW signal is 
"mounted" over a larger scale SW mode (of negative ampli- 
tude) it can produce a net negative contribution to wtg- 

Our findings are in agreement with recent work on the 
cross-correlation measure of WMAP with NRAO VLA Sky 
Survey radio source catalogue (NVSS) (Boughn & Critten- 
den 2004), (Nolta et al 2003). They detect a signal of 
wtg — 0.16/aK with 1.8 deg pixels (148 counts/pixel), what 
is consistent with our measurements once the different se- 
lection function is taken into account. Our cross-correlation 
analysis on large-scales reveals that the evolution of the 
gravitational potential has been strongly suppressed at low- 
z. The drop of this positive correlation on small-scales sug- 
gests that we might be measuring the SZ-induced distortion 
of CMB photons by nearby clusters. Deep large area galaxy 
surveys, such as the SDSS, should be able to confirm these 
results, provide tighter constraints on cosmological param- 
eters and improve our knowledge of cluster physics (Peiris 
& Spergel 2000). Such analysis, together with a detailed 
treatment of the SZ and lensing effects will be presented 
elsewhere (Fosalba, Gaztanaga & Castander 2003). 
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